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Among the 0{as tan/3) contributions to tiie Wilson coefficients C? and Cg, relevant for the decay 
B —> Xs7, those induced by two-loop diagrams with charged-Higgs-boson exchange and squark- 
gluino corrections are calculated in supersymmetric models at large tan/?. The calculation of the 
corresponding Feynman integrals is exact, unlike in previous studies that are valid when the typical 
supersymmetric scale Msusy is sufficiently larger than the electroweak scale mweak(~ mw,mt) 
and the mass of the charged Higgs boson mj^± . Therefore, the resufts presented here can be 
used for any vafue of the various supersymmetric masses. These resufts are compared with those 
of an approximate cafcufation, afready existing in the fiterature, that is at the zeroth order in the 
expansion parameter (n-^gj^j., m|f )/M|usy) s^nd with the resufts of two new approximate cafcufations 
in which the first and second order fn the same expansion parameter are retained, respectivefy. This 
comparison affows us to assess whether the resufts of these three approximate cafcufations can be 
extended beyond the range of vafidity for which they were derived, i.e., whether they can be used 
for niH ^ Msusy and/or Msusy ~ rn^aa,k- It is found that the zeroth-order approximation works 
weff even for tuh ^ Msusy, provided M|ugY ^ 'Tiwcak- The incfusion of the higher-order terms 
improves the zeroth-order approximation for m'jj ^ M|usy) but it worsens it for tuh ^ Msusy- 

PACS numbers: 12.60. Jv,13.20.Hc,14.80.Cp 



I. INTRODUCTION 



While the rate for the inclusive decay B — > Xs^, BR(S Xsj), has been calculated up to the next-to- leading order 
(NLO) in QCD within the standard model (SM) ^Ij, similarly precise calculations exist for only some extensions of the 
SM. Achieving a NLO precision for this decay rate in such extensions, and in particular in supersymmetric models, 
is important. It seems unlikely that an increased experimental precision in the three apparatuses where this decay 
is measured, i.e. BABAR, BELLE and CLEO 0, combined with the increased theoretical precision of a possible 
estimate of the next-to-NLO corrections in the SM |^, will bring unequivocal signals of new physics. Nevertheless, 
calculations of BR(5 X^j) with NLO accuracy for the models that are considered the most likely candidates to 
extend the SM could help in understanding where the effective scale of these models sets in, and the extent of the 
spreading of masses of additional particles around this scale. 

NLO calculations exist for two-Higgs-doublet models of type I and type II ^3, 'sl, '6l| , as well as models in which 
the couplings of the charged Higgs boson to fermions are, in absolute values, those of type II models, but are, in 
general, complex 0]. In the case of supersymmetric models, the situation is as follows. For generic models, the QCD 
corrections to the electroweak rate, calculated first in ?], have been included only at leading-order (LO) precision j^. 
Higher-order QCD corrections have been evaluated for specific scenarios. 

In one class of such scenarios 9j, the two charginos, one t squark, which is predominantly right handed, and 
the charged Higgs boson are assumed to be relatively light, while all the other squarks and the gluino are heavy. 
Moreover, no additional sources of fiavor violation are present at the electroweak scale, other than the Cabibbo- 
Kobayashi-Maskawa mixing elements. 

Other studies [ToL[lTI| considered a slightly different class of scenarios, in which tan/3 is large, and the supersymmetric 
spectrum is like the spectrum for the scenarios of Ref. |^ , but without the assumption of a light i squark. The same 
minimality in fiavor violation is also assumed |12| . The importance of supersymmetric models at large tan (3 hardly 
needs to be highlighted here, as tan/J tends to be large whenever an attempt to unify Yukawa couplings is made, 
as required by a SO(IO) grand unification Indeed, many phenomenological calculations exists, in which the 

constraints imposed by the measured rate of i? ^ Xg'j to models at large tan/3 have been analyzed at LO in 
QCD The papers in Refs. [Taiia allow one to refine these analyses in models which, as well as having large 
tan /?, predict the same type of supersymmetric spectrum that they assume. 

The restriction to specific ranges of masses of the supersymmetric partners of the SM particles, made in Refs. 0, 
[Tol , has made possible the use of an effective Lagrangian formalism, and has led to rather compact and simple 
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formulas. These very same restrictions, however, hmit the usefulness of such calculations, unless it is proven that the 
formulas obtained in these papers can be safely used beyond the range of validity for which they have been derived. 
Nevertheless, they have often been employed also when the charged Higgs boson is as heavy as the squarks and the 
gluino. It is interesting to understand, in such cases, how far the resulting analyses are from ideal ones, in which 
exact formulas are used. 

In this paper we address such an issue, focusing on the scenarios of Refs. [lol[Tl| . We consider the gluino-induced 
supersymmetric corrections of 0{as tan/3) that are the largest beyond-leading-order corrections in scenarios with large 
tan (3. Corrections of this type are obtained by (i) including the finite corrections to the 6-quark mass in the fermion- 
fermion-Higgs-boson and in the fermion-sfermion-Higgsino couplings |l5l IT^ , and (ii) "dressing" with squark-gluino 
subloops the one-loop diagrams mediated by the charged Higgs and the charged Goldstone bosons that contribute to 
B — > Xsj at the partonic level. A graphical representation of this "dressing" is explicitly shown in Fig. ^ for the 
diagrams contributing to 6 — > S7 or to 6 sg and mediated by the charged Higgs boson. (In this figure, the photon 
or the gluon is assumed to be attached in all possible ways.) Indeed, the substitution of the vertex SLtRH^ of the 
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FIG. 1: 6 — > sy/b — > sg by charged-Higgs-boson exchange, at one loop, on the left side, and with inclusion of gluino-squark 
subloops, on the right side. The photon/gluon is assumed to be attached in all possible ways in both diagrams. 

one-loop diagrams with s^IrH^ , possible in the two-loop diagrams, allows the gain of a tan/3 factor, yielding the 
required as tan /?. 

As already pointed out, in the calculations of Refs. [Tol ITU, all scalar superpartners of the SM fermions, at the 
scale MgusY, are assumed to be heavier than the t quark and the W boson, whereas the charged Higgs boson is 
assumed to be around the electroweak scale mweak, with mweak ~ Mw, rnt- In the limit of large AfgusY, the two- loop 
diagrams in Fig. ^ in which the photon is emitted only by the t quark and the charged Higgs boson, and the gluon 
only by the t quark, all with chirality flip on the t-quark line, are of nondecoupling nature. In the same limit, all other 
diagrams in Fig.Qldecouple. It is conceivable, therefore, that for r7i^± ~ to^ ^ -^susY' ^bc former two-loop diagrams 
are the only ones that give a sizable contribution to the 0{as tan (3) corrections to the decay amplitudes of 6 ^ 57 and 
b — > sg, and that the expression for their squark-gluino subloops is well approximated by the expression for the same 
subdiagrams with vanishing external momenta. Thus, in this approximation, these two-loop diagrams are factorized 
into two one-loop diagrams. Such a factorization is supported by the use of an effective two-Higgs-doublet Lag rang ian 
formalism, in which all heavy degrees of freedom are integrated out. (A similar formalism was used in Ref. |l7|. in 
which corrections of 0{a tan f3) to the decay amplitude of 6 — > 57 were calculated.) 

A potential problem with this approach may arise when the supersymmetric spectrum is not much heavier than 
Wwcak, and/or when the mass of the charged Higgs boson tends to be closer to MgusY than to Wweak- (Papers with 
further improvements of the original calculations, such as those in Ref. ISj, do not address this issue.) One plausible 
solution to this problem is to include in the original calculation higher-order terms, up to 0{{m1,^^^^, ™h± /-^susy)")' 
with a suitable n. That is equivalent to saying that the effective Lagrangian is extended to include higher-order 
operators. (A discussion on this point can be found in Ref. |0|.) An efficient way to carry out this extension, 
consistently including all the operators that yield terms of the same order n, is to make use of the heavy mass 
expansion (HME) 20J. Starting, for example, from the regime in which m^,M^^ ^ ^ -^^susyi ^^i^ technique 

allows one to add all operators needed when mjy± tends to MgusY from below. It remains, however, to be established 
up to what value of n it is necessary to extend this expansion, in order to obtain an estimate of BR(i3 Xgj) 
adequate for all possible values of Msusy and of predicted by different supersymmetric models at large tan (3. 

Clearly, an accurate answer to this question can be given only by a comparison with the exact calculation of all the 
two-loop diagrams that give rise to 0{as tan (3) corrections, in which no assumption is made on the relative size of 
TOjj±, Msusy, and mweak- 

In this paper we perform such an exact calculation of these two-loop diagrams contributing to the two Wilson 
coefficients Cj and Cg related to the partonic processes b ^ sj and b sg. Among all the contributions to these 
processes we restrict ourselves to those mediated by the charged Higgs boson, addressing first the problem of 
obtaining results valid throughout all values of m^± . We postpone the presentation of the complete BK{B Xg'j) 
to future work (2 a] . 
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We also calculate the same diagrams using the HME technique, up to order 0{{ni^^^^,m'j^j^ /M^jj^y)"^). The hope 
is that, in the case in which the series expansion appears to be quickly converging, approximate and possibly still 
compact formulas can be provided, which nevertheless allow one to extend the results existing in the literature to all 
possible values of to^± and Msusy- It turns out that the series expansion in to^± /Afgugy) is far from being 

as well behaved as hoped and that, in general, the use of the exact results cannot be avoided. However, whereas the 
approximate calculations including terms up to 0((m^^^j,, M'^susy)^) S^^^ results in disagreement with the exact 
one for values of mjy± ^ Msusy, the lowest-order term of the series expansion, i.e., the nondecoupling approximation 
of Refs. ^ili, seems to be a good approximation of the exact result throughout all ranges of m^± , provided AfsuSY 
is sufficiently larger than mwcak- Deviations appear for supersymmetric particles not much above the electroweak 
scale. 

The paper is organized as follows. In Sec. ^ we list all the diagrams needed to evaluate the charged-Higgs- 
boson mediated 0(as tan /3) contributions to the Wilson coefficients Cj and Cg entering in the calculation of the 
BR{B Xg'^). We also give these coefficients at the matching electroweak scale fxw in terms of two-loop scalar 
integrals. In Sec. Illll we review some basic properties of the HME technique and discuss how it can be applied to our 
specific case. In Sec. lIVI we show numerical results for the charged-Higgs-boson contributions to C7{fiw) and Cs{fJ,w), 
obtained from the exact two-loop calculation, and from the two approximate calculations including up to the first two 
orders in {'m'^^^^^,m'jj±)/Mg^gY. These results are compared and discussed. Finally, in Sec.|3we conclude. 



II. THE TWO-LOOP DIAGRAMS 



The charged-Higgs-boson mediated two-loop diagrams, contributing at 0{as tan/3) to the partonic decays b ^ sj 
and b — > sg, are obtained from the right diagram of Fig. Rafter allowing a photon and/or a gluon to be emitted in 
all possible ways. They are shown explicitly in Fig. 12 where the photon may be replaced by a gluon, and vice versa, 
whenever possible. 

It should be observed that, while the one-loop diagram on the left side of Fig. ^ has a chirality flip in the internal 
t-quark line (we neglect in this discussion the diagram with chirality flip on the external fo-quark line, which is tan /3 
suppressed), the diagrams in Fig.|21can have such a chirality flip on the t-squark line also. The number of contributions 
to be evaluated, therefore, amounts to 8 for the calculation of the Wilson coefficient C7{nw), and 8 for the coefficient 
CsCmw)- The diagrams used for the calculations of BR(_B — > Xg-f) in Refs. JlJ are the two at the top of Fig. [21 
for the calculation of Crdiw), i.e., Ufa) and[2Ib), and the diagram|2fa), with the photon replaced by a gluon, for the 
calculation of CsifJ-w), all with chirality flip on the t-quark line only. 

We have calculated all the 16 contributing terms, by making use of results and techniques presented in Ref. 0|. Our 
normalization of the Wilson coefficients Ct(iiw) and CsifJ-w) is the conventional one, as follows from the definition 
of the effective Hamiltonian, 

TJeff D --^V,:Vtt [C7(m)07(m) + Csifi)Ositl)] , (1) 

and of the operators O7 and Og, 



where and GJ^^ are the field strengths of the photon and the gluon, respectively. We denote by Cr^nifJ'w) and 
Cg.Hip-w) the tan /3-unsuppressed charged-Higgs-boson contribution to Cj^fiw) and Cs{fiw)i and we decompose them 
as 

C^,H{^^W) = . ^ , V ^ [C°hM + AClHit^w)] , (3) 
1 + Af,^,6tan^ 

where C^j^diw) and AC} j^{fiw) are induced by the one-loop diagram in Fig. and the two-loop diagrams in Fig.|21 
respectively. The overall factor 1/(1 + Ab^.btan/3) (see notation of Ref. ^^) stems from expressing the H^t^bji 
Yukawa coupling in terms of rrib, corrected up to 0{as tan/3) jl5i I1(t||. Atj^^b is given by 



where the function / is 
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FIG. 2: Charged-Higgs-boson mediated diagrams contributing at order 0{as tan/3) to the partonic decays & — > 57 and b — + sg. 
The photon must be replaced by a gluon, and vice versa, whenever possible. 



2 2 2x mfml ln(mf /to^) + m^m| ln(m^/m|) + TO^mf ln(TO§/TOj) 
As the photon can be emitted by the t quark and the charged Higgs boson, we have two contributions to C^{^\y)- 



while only one to Cg{^w)' 



C'^hM = \Qh^fJ^, (6) 



ct:hM = A^F, (-^) , (7) 



2 V"ff±, 



due to the emission of the gluon from the t quark. The superscript indices a and h in the above expressions denote 
emission of the photon/gluon from the t quark and the charged Higgs boson, respectively. The functions and 
are listed in Appendix^ Qt and Qh indicate the electric charges of the t quark and the charged Higgs boson 
{Qh = -1). 

The contributions to AC^ jLf(/iiv), due to the different two-loop diagrams in Fig.|21 are 

AC]%{ixw) = ^QtCp {as tan/3) mt {4nf [{Ut)2^{U^)^2 mtm~g Iu2 - {Ut)2,{Uj')^iIui] , 

^C^^i^iw) = ^ QhCf (as tan/3) mt fi {Anf [{Ut)2z{U^)t2 mtm~g I„i2 - {Ut)2iiU^)tiImi] , 

ACl'^ifiw) = I Qt Cf («. tan/3) mt ^ [Anf [{Ut)2^{UJ)a mtm-g If,^ - {Ut)2^{U^)^lIul] , 

ACljjM = \ Qs Cf {as tan/3) {Anf [{Ut)2^{U^)^2 rutmg ha - {Ut)2^{UJ)^lIMl] ■ (8) 
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In the above expressions, the definitions 



tan/3 = ,2^%+^ (9) 

VD 2 



are adopted, with vu and vd the vacuum expectation values of the neutral components of the Higgs doublet with 
hypercharge +1/2 and —1/2. All phases for the supersymmetric parameters are assumed to be vanishing. Thus, the 
matrix Ut is the 2x2 orthogonal diagonalization matrix of the matrix for the t-squark mass squared, and a summation 
over the index j, identifying the two t eigenvalues, ti and t2, is understood. Our conventions for the squark sectors are 
listed in Appendix^ We have also assumed that the left-right mixing in the matrix for the s-squark mass squared 
is nearly vanishing and therefore only one s-squark eigenstate, the left-handed one, denoted simply by s, is needed 
for this calculation. Qi and Qs indicate the electric charges of the i and the s squarks. Cp = 4/3 is a color factor. 
Finally, the scalar integrals lui, Iti2, Ihh, lHi2, lui: hi2^ -^sn, and Isi2 are listed in integral form in Appendix [HI 

Similarly, the contributions to ACg^(/ivi') are 

ACl'^Hif^w) = ^Cf {a,tsinp)mt fi {4TTf [{Ut)2t{U^)t2mtm~g Iti2 - {Ut)2t{UryaItii] , 
AC^X/iv^) ^ ^Cf (l-^) (as tan/3) TO* A' {4nf [{Ut)2;{U^),2mtmg 1^^, - {Ut)2^{Ul)^lIi^l] , 
ACi^(m) ^ (l-^) (astan/3)mtAi i^^rf [{Ut)2^{Uf),2mtmg ha - {Ut)2^{Ul)^lh^l] , 

^Cl^nim) ^ \Cv {as tan/3) ^l (4^)^ [{Ut)2^{UJ)^2 mtuig Ig,2 - {Ut)2^{Uf)^lIg^l] , (10) 

where, obviously, AC| ^ (/iyt/) — 0. Here, Cv = 3 is another SU(3) group factor and the integrals Igii,Igi2 are also 
defined in Appendix O The explicit expression of the integrals introduced in this section can be obtained upon 
request as a FORTRAN code. 

The results shown in Refs. [lillll for the BR{B X^j) are obtained using only the first term in the square brack- 
ets of ACj°^{iiw), If (/xi4/), and ACg|f(/in/), with the approximations for Iti2 and lHi2 listed in Appendix lUl 
As anticipated in the Introduction, we call the approximation of these references the nondecoupling approxima- 
tion. It is indeed of nondecoupling type in the limit of heavy supersymmetric particles and collects all terms of 
0((m^g^j., m^i/MgugY)"), as explained in greater details in the next section. Strictly speaking, however, it contains 
formally decoupling terms coming from the masses and the couplings of squarks (see Ref. llSj). 



III. EFFECTIVE LAGRANGIAN AND HEAVY MASS EXPANSION 



The calculation of the two-loop diagrams discussed in the previous section was performed in Refs. |lOl lllj under 
the assumption of heavy squarks and gluino, at the scale Msusy, and a light charged Higgs boson, of O(mwcak) 
(wweak ~ Mw , mt) , with ml^^^^ ^ -^^susy- III these calculations, terms of 0((m^^^j., m^i/Mgugy)") with n > 0, 
i.e., terms that decouple in the limit Mrusy — > c», were neglected except for those trivially accounted for by the 
couplings of squarks and their masses [l8| . 

The amplitude of the diagrams in Fig. |21can be expanded in (w^^j^^]^, to^± /Mgygy). It seems plausible that by 
retaining in such an expansion terms up to ^((m^^jjjj., to^± /Afgugy)") with n > 1, may allow one to extend the 
results obtained to values of toh± not too dissimilar from MgusY, as well as to a not too large value for AfgusY- In 
this spirit, it seems worth performing these approximate calculations, in the hope that the results turn out to be still 
relatively compact. 

A systematic way to include all the needed terms at 0((m^j,^j^, m^i/Mgugy)") is provided by the heavy mass 
expansion technique pol |. Before proceeding further, we recall here some basic properties of this expansion. 

Let us assume that all masses of a given Feynman diagram F can be divided into a set of large M = {Mi, M2, . . .} 
and small m = {mi, m2, . . .} masses. If all external momenta q = {qi, 52, ■ • ■} are small compared to the scale of 
the large masses M, then the dimensionally regularized (unrenormalized) Feynman integral Fr associated with the 
Feynman diagram F can be decomposed as 

Fr ° '^g-',rn-'F^ {q\m\M) , (11) 
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FIG. 3: The diagrams r/71, on the left side, obtained from the diagram (b) in Fig. |21 The corresponding 71 are shown on the 
right side. [See Eq. IIH in the text.] The chirality fiip is on the t-quark line in the upper two diagrams, on the f-squark line on 
the lower two. All diagrams in this figure and in Figs. 0] and |S] below are shown in the gauge eigenbasis for squarks and Higgs 
bosons. 



where the sum is performed over all subdiagrams 7 of F which (i) contain all lines with heavy masses (M) , and (ii) 
are one-particle irreducible with respect to lines with small masses {m). The case 7 = F is always included in the 
sum For each 7, rrf' denotes the set of light masses, q'' the set of all external momenta with respect to the 

subdiagram 7, which can be internal momenta with respect to the full diagram F. The operator T performs a Taylor 
expansion in the variables q"^ and rnP' and it is understood to act directly on the integrand of the subdiagram 7. The 
diagram F/7 is obtained by reducing 7 to a vertex in F. Thus, by factorizing the product of scalar propagators of the 
original, say, /-loop diagram F as Ilr ~ Xlr/^II^, the decomposition of the original Feynman integral Fr is simply 

Fr/y oTqi^m-rFj = J dki ■ ■ ■ dkl Ur /-y Xp ^m-r^-i ■ (12) 

Note that the Taylor operator T introduces additional spurious IR or UV divergences in the various terms of the 
sum which cancel in the sum. 

For the calculation of the diagrams in Fig. |21 we start assuming that 

"^ff± ~ m^oak < ^^SUSY, (13) 

where mweak is, in turn, 3> ruh. In this limit, the diagram (b) in Fig. [3 with chirality flip on the i-quark line has a 
Feynman integral that can be decomposed into two terms contributing to the sum J^-y of Eq. In the first term, 

71 is the upper right diagram of Fig. 13 and the corresponding Feynman integral must be expanded in the momenta 
of the external particles, of the charged Higgs boson H^, and of the t and s quarks; r/71 is the diagram on the 
upper left of the same figure. Although the calculation is actually done in the mass eigenbasis for squarks and Higgs 
bosons, the diagrams in this figure (as well as those in Figs. 0] and [S]) are shown in the gauge eigenbasis. In an effective 
two-Higgs-doublet Lagrangian, obtained after integrating out all heavy degrees of freedom, i.e., squarks and gluino, 
the zeroth-order expansion of 71 is described by the term in the Lagrangian 

Z:^^ = -V.s'^^^A,,, H+tnSL + (H.c), (14) 
V sm fj 

where, again, we follow the notation of Ref. ^^'^ coefficient Aj^^s is 

C (X 

At„,s = -^--^fim~g{Ut)2z{Ui^)z2 I{ml,mj,ml) . (15) 

The function / is the same as in Eq. This zeroth order corresponds to the nondecoupling approximation of 

Refs. and gives rise to the following contribution to the Wilson coefficient Cr^HifJ-w)- 

1 1 TTh^ [ TTXt' 

^^•7 (/iw/)|„ondcc = - -QifCi^ (as tan /5) —(\Jt)2i{Vj)i2 I{m\,m\,'rr?g) J Fi[—-^ 

z Zn » -''^H± V^//±, 
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FIG. 5: The diagrams r/71 for the transition b sg, on the left side, with the gluon emitted by the gluino g. The corresponding 
71 are shown on the right side. The chirahty flip is on the t-quark line in the upper two diagrams and on the t-squark hne on 
the lower two. 



The second term contributing to the sum "^^^ 72, is the full diagram, and its Feynman integral is to be expanded 
in mt,TO^± jTOb, and the external momenta; r/72 is, trivially, the identity. Again, in an effective two-Higgs-doublet 
Lagrangian, the operator describing 72 is the very operator O7 in Eq. (0). This gives rise to a vanishing zeroth-order 
term in the expansion o f lig ht masses, i.e., a vanishing nondecoupling contribution, and was, therefore, not included 
in the analyses of Refs. [iflllj- 

The case in which the chirality flip is on the i-squark line instead of on the i-quark line can be treated in a similar 
way. In this case, the term 7^ shown by the lower right diagram of Fig.|31 is already of decoupling type and was also 
not considered in Refs. [lOlllll|. 

The procedure to be followed is similar in the case in which the photon is attached to the t-quark line. When the 
chirality flip is on the t-quark line, the diagram 71 is again the upper right diagram of Fig.|3| which at the zeroth-order 
expansion gives rise to the same term of the effective two-Higgs-doublet Lagrangian given in Eq. H14(l . The diagram 
r/71 is analogous to the upper left diagram of Fig. |21 The corresponding contributions to the Wilson coefficient 
Cji^iw) and to CsilJ-w), when the photon is substituted by a gluon, are 

1 TTl^ ( TtT/^ \ 

ACl'k{^iw)\non<iec ^ ^Qt tan/3 A*,,, , (17) 

Again, the contributions from the second diagram 72 and the diagrams obtained when the chirality flip is on the 
t-squark line are of decoupling type. By comparing Eqs. (|16() . 117() . and 1181) with the tan /3-unsuppressed one-loop 
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contributions jj{^w) in Eqs. © and lO, the relation 

^CIh {p.w)\nondcc = - tan/3 Af^^^ C^,H{^^w) (19) 
follows, which shows clearly that /S.CI fj{iiw) has the same mt and rafj± dependence as that of fj{^w)- 

Slightly different is the case in which the photon is emitted by one of the two squarks t,s. In the case of emission 
from the s squark, for example, the decomposition of the Feynman integral has, again, a term corresponding to 
the full diagram F^^, to be expanded in mt,mH± jmb-, and the external momenta. The other term is given by the 
convolution of the Feynman integral corresponding to the upper or lower box diagrams 71 shown on the right of 
Fig- El (depending on where the chirality is flipped, i.e., on the t quark or i squark), expanded in the momenta of H^, 
of the t and the s quarks and the momentum of the photon, and the Feynman integral -^1^/71 of the upper or lower 
diagrams r/71 on the left side of the same figure. 

When it is a gluon to be radiated off during the interactions that lead the b quark into an s quark, there is 
the additional possibility of having the gluon emitted from the g line. The decompositions in 71 and r/71 of the 
corresponding diagrams are explicitly shown in Fig. |S| 

The results of the calculation in the regime of Eq. H13|l turn out to be rather involved, already at 
'^('^wcak' "^ffi/^susv) ^'^d 0((m^^jj]^, m^i/Mlugy)^)- We report them only numerically in the next section, to- 
gether with the numerical results of the exact calculation of Sec. |n] 



IV. H+ CONTRIBUTION TO Cr^ifJ-w) UP TO C>(q, tan/3) 



We present here numerical results for the charged-Higgs-boson contributions to the Wilson coefficients C7 and 
Cg at the scale /x^/ = M\y, i.e., C'!^h{iJ'w) and C%^h{.I^w), including the 0(as tan/3) corrections discussed in the 
previous sections. We make a comparison of the exact results for Ct^hIpw) and CsMifJ-w), which we denote by 
C7.H{^J■w)\cxl^ct and C8,_f/(/^w)|cxact, obtained from Eqs. (jHJ and H1(J|) . with several approximate forms, C7^h(/^w)| approx 
and Cg^ij (/ivF)|appiox: the nondecoupling approximation, and the two approximations in which the first- and second- 
order terms in the expansion of the exact coefficients in (jn^^^^, "^lf± )/-^SUSY are retained. 

In Figs, ini and [7| we plot the ratios 

/ . _ Ci^nifJ'Wjlapprox — Ci^H{f^w)\ciLact ,. 

rt(pw)^ 7 r. (1 = 7,8), (2U) 

^i,H[f^W Jlexact 

in which the 0{as tan/3) corrections to the 6-quark mass cancel out, showing the goodness of each approximation as 
a function of toh± • We denote these ratios for the nondecoupling approximation and those with an expansion in the 
first and second order in (m^ , )/M|ugY by dotted, dashed, and dot-dashed lines, respectively. 

Two sets of parameters for supersymmetric particles are used. For Fig.El we have chosen a superpartner spectrum, 
called here spectrum I, with {ms^,mQ,mfa,mga) = (700, 450, 435, 470) GeV, At = 150GeV, mg — 600GeV, and 
[X = 550 GeV. For Fig. [7| a lighter spectrum is considered: (mg^ , Wg , ra j,,, , m-^c ) = (250, 230, 210, 260) GeV, At = 
70 GeV, nig — 200 GeV, and /i — 250 GeV. This is denoted as spectrum II. As for other input parameters, we have 
used tan/3 = 30, mt{^J.w) = 176.5 GeV, which corresponds to a pole mass Mt — 175 GeV, ■mb{^w) = 3 GeV, Ab = 0, 
a^inw) = 0.12, Mz = 91.2 GeV, and sin^ Ow = 0.23. 

For spectrum I, the difference between the exact calculation and the nondecoupling approximation is very small 
in the whole range of to//±, even for rn//± ^ Msusy- This is an unexpected result since, as discussed in Sec. IIIII 
the nondecoupling approximation is, in principle, theoretically justified only for ?7i^± ^ -^susy- 1^- the case of 
spectrum II, ry^g become larger for the nondecoupling approximation. The corrections beyond this approximation are 
of the same order as the SU(2)xU(l) breaking effects in the supersymmetric particle subloops and are no longer 
negligible. Nevertheless, ry and rg for the nondecoupling approximation remain of the same order of magnitude for 
increasing mu±, up to mjLf± ^ Msusy- 

In both cases, the bulk of the difference between the results of the exact and nondecoupling calculations come from 
Iti2 in ACy 5/ and ACg 5/, and also from Igi2 in ACg of Eqs. © and (fTUIl . The inclusion of the (rn^ , A^susy)" 
terms by the HME, as described in Sec. IIIII improves the goodness of the nondecoupling approximation when 
TO//± < AfsuSY, but worsens it when TO^f± <^ .^^susy- This is clearly seen in Fig.Q 

In Fig. |S1 we show C^^H(^J^w) and Cf,^H{l^w) for the nondecoupling approximation and exact two-loop results, as 
well as the one-loop results C'°^(/x^), for spectrum II. The 0(as tan /3) corrections are comparable to the one-loop 
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FIG. 6: Ratios r7^s{l-i-w) for various approximations of the 0(astan/3) corrections, as defined in the text, for spectrum I. The 
dotted, dashed, and dot-dashed lines show the nondecoupling approximation, first-order expansion, and second-order expansion 
in (m?,m^±)/M|usY, respectively. 
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FIG. 7: r7^s{^^w) for spectrum II. Notation is the same as in Fig.|S| 



contributions. The deviation of the exact calculation from the nondecoupling approximation is small, but nonnegligible 
for Cs,H- It is also shown that all three results have a similar dependence on mjj±. 

To understand the results for mjj± ^ AfgusY qualitatively, we focus on the diagram in Fig.jS^a), with chirality flip 
on the t-quark line. When mfj± is sufficiently larger than mt, the contribution of this diagram to ACj jj{^w) f^nd 
AC^ jj{^w) is among the largest. Analytically, it is proportional to iJ,mgIti2, with the integral Iti2 listed in Eq. (|C3|I 
of Appendix [ni For the following discussion it is convenient to write fj.'mglti2 in the form 



fJ-mgIti2{mt,mH± , mf.,ms, rUg) 



(27r)4 [P„to2]3 [fc2_^2^^] 



Yu2 {P;m^,,ms,mg) , 



where Yu2{k'^',mj:. ,ms,mg) represents the subdiagram contribution to the vertex H SLtR and is given by 
Yu2ik'^;mf.,mg,mg) = unig [-2F + {k^ - m1)G\ (jt'^;m?^,rr>,j,mfj , 



(21) 



(22) 
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FIG. 8: C7,£r(/ivi/) and Cs,h{ij-w) for spectrum II. The dotted, dashed, and sohd hnes show the one-loop result, nondecoupling 
approximation, and exact two- loop result, respectively. 



with 



CI// 2 2 2 2\ 

F{k ;mr,TOj,ms) 



dH 






1 




(2^)4 








P — TO? 


dH 






If' 




(2^)4 




- fc)2 — m? 




;2 _ jy^2 



2 • 



(23) 
(24) 



The nondecoupling approximation /iTOg/ii2 Inondoc of Eq. (|21|l . where the integral /ti2|noiidoc is given in Eq. (|C13|I . is 
obtained by approximating Yti2{k'^', fni., rris, nig) by 



(25) 



which is a constant with respect to k^. Note the relation FlQ-jni- ,m'i,m'i) = (— i/167r2)/(m? , to|, to?), where / is 
the function defined in Eq. Q. To simplify our discussion of the behavior of these functions, we hereafter set mj. , 
rris, nig, and /i equal to Msusy- 

For < M^ygy, F{k^;M^^gy) and G{k^; M^jjgy) behave as 



F{k';Mi^SY) 



O 



G(/c2;M|usy) = O 



^"SUSY 



-'"SUSY 



O 



M4 

^'^SUSY 



(26) 



For |fc2| -^'j^susY' the behavior is 



Fik';Mi^sy) ^ 

G(fc^;M|ugY) ~^ 
The behavior of Yti2{k^; M^^^y) is therefore 

^ti2|nondcc + O 



k'' M, 



O 



1 



In- 



SUSY 



-^SUSY 



(27) 



i^t.2(fc'; Msusy) 



fc2 



to; 



^'-'SUSY -""SUSY 

O ( Mjsy 1^ fc' 



fc2 



Af 2 

-""SUSY 



(|fc'l«M|usY), 
(|fc2|»M|usY), 



(28) 



11 



which supports the naive expectation that a substantial deviation of Iti2(jnt,mu± , M'^^o^y) from 
/ti2(™t,m-ff±,^|usY)lnondec may arise from the region \P\ 'k, M^^^y. 

The factor multiplying Yti2{k'^; -^susy) ^ Eqs. (|21|l . however, plays a rather important role, leading to the fact that 
this expectation does not hold in the case in which A/susy is not rather light. Since for |fc^| 3> rrijj± this factor drops 
as d^k/k^, the integral 1)2111 gets its largest contribution from the region |fc^| ^ m?^±. A closer inspection actually 
shows that it is the region of small up to |fc^| — 0{mf), which determines the bulk of the value of this integral. 
If rrit is sufficiently smaller than Msusy, ^ti2(fc^; -^susy) "^o^^ ^'^^ deviate substantially from lti2|nondcc- The region 
in fc^ where the deviation Yti2{k^ ; Mg^^y) — yti2|nondcc is largest, i.e., |fc^| ~ -^^susY' weighted by a rather efficient 
suppression factor in Eq. (|^ if mf <^ ^ -^Iusy- 

It is clear that an expansion of the integral Iti2{mt, mfj± , Afgugy) in "m-t / M^u^y and m^± /-^susyj obtained from 
the HME, generates terms 0{mP^± ) at the first order, and terms 0{'m?^± ) at the second order. These terms contribute 
to give a better approximation of the exact results for Ci_//(/in/), when to^± < Afgygy, but have a dependence on mfj± 
that is rather different from that of the exact results, when to^± ^ ^Iusy- already observed, the m}j± dependence 
of the exact results for Ci^uilJ'w) a-t order 0(as tan /3) is similar to that of their nondecoupling approximation, which 
is, in turn, identical to that of the one- loop results C'^^{fj,w) (i-e., the results at LO in QCD) for the same Wilson 
coefficients. As for the series expansion obtained through the HME, the qualitative discussion sketched above and 
the limited number of terms we have calculated do not allow us to conclude whether it is convergent or not. Even if 
convergent, however, the numerical results we have plotted in the previous figures show clearly that it is not converging 
fast enough to be of any practical use in the region mj^^ ^ -^susY' originally hoped. 



V. CONCLUSION 



In this paper we have evaluated the 0{as tan/3) corrections to the charged-Higgs-boson mediated contributions to 
the Wilson coefficients relevant for the decay B Xgj, in supersymmetric models with large tan/3. These corrections 
are generated by the shift of the 6-quark mass in the Higgs-boson-quark couplings and by the dressing of the one-loop 
charged-Higgs-boson diagrams with squark-gluino subloops (see Fig.|2l. The former corrections are very well known. 
In this paper we have focused on the latter class of corrections, which exist in the literature in an approximate form. 

In previous studies 0, ^| , the contributions from these two-loop diagrams were calculated using an effective two- 
Higgs-doublet Lagrangian formalism where the squarks and gluino are integrated out. This method is theoretically 
justified in the limit in which the charged Higgs boson is light, i.e., at the electroweak scale ~ Wweak, whereas the squark 
and gluino masses, ~ AfsuSY, are rather large. The resulting approximation, in which only all the nondecoupling 
terms in the large MgusY limit are retained, has clearly the nontrivial advantage of providing a rather compact 
result for these corrections. However, it is in general expected that its deviation from the exact two-loop result, of 
<^{'^weak''^ff±/-^SUSY)' bccomcs significant for Wwoak ~ Msvsv and/or m^± > Msusy- 

We have calculated the contributions of the two-loop diagrams in Fig. Inexactly, without assuming any patterns for 
the masses of the particle involved. By making use of the heavy mass expansion technique, we have also evaluated two 
additional approximate forms for the same diagrams, including all terms up to 0(m^g3^jj, r7i^± /Mlygy) in one case, 
and all terms up to 0{{m'^^^^^,m'^±/Mg^gy)^) in the other. The exact calculation, not confined to specific values of 
the masses involved, has allowed us to establish the goodness of all three approximations. 

Our findings can be summarized as follows. Surprisingly, the results of Refs. [Tollll| | approximate the exact two- loop 
result quite adequately, irrespective of the value of mH± relative to MgusY, provided Msusy is large enough. The 
unexpected absence of a large deviation for the case of m^± ^ AfgusY with rn^^^^ ^ -^susy understood from 

the structure of the two- loop integrals. This points to the fact that the only relevant hierarchy in this problem is that 
between Wwoak and Msusy, whereas, for m^± > rriwcak, the value of with respect to that of Msusy plays a 

rather marginal role. Therefore, deviations between the exact result and that of the nondecoupling approximation can 
be found only for Msusy not much larger than mwcak- The inclusion of the high er-order terms 0(rn^^^^, /-^susy) 
and 0((to^j,3^]^, TO^i/Mgugy)^) improves the approximation of Refs. 0, 0| for m^± <C Mgygy, as expected, but 
makes it worse for to//± <; Msusy- This behavior is attributed to the fact that these higher-order terms in the 
(™wcak' "^ffi M'^susy) expansion have a dependence on m/^± rather different from that of the lowest-order terms in 
this expansion (the terms of the nondecoupling approximation) and that of the exact two-loop result. 

We have illustrated our findings by showing the values of the Wilson coefficients C7 and Cs at the electroweak 
matching scale /iw, for different gluino, squark, and charged Higgs boson spectra. We postpone a presentation of the 
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same coefficients C7 and Cs at a low scale ~ irib and of the actual branching ratio BR(i3 ^sj) to future work [23 ■ 
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APPENDIX A 

We report here the functions ^3(2;) and F4{x) introduced in the second paper of Ref. [2]| : 

Fi{x) = ^^^^^(a;2-l-2xlogx). (Al) 



APPENDIX B 



We have adopted the following conventions for squark mass eigenstates qi ,q2 and mass eigenvalues m?^ ,m?^ : 

where the diagonalization matrix Uq is such that 

U^MlUq^Ml^^< . (B2) 

Al| is the mass squared matrix for the squark q and Al~ is the diagonalized squark mass squared. For q — b and i, 
the matrix | is, respectively, 

^ ( mb{Ab- n tan (3) \ 



and 



with Dj^ ^ and D)^ ^ given by 



A^2^( ™| + "^t+^i rnt{At- cot (3) 

* ' nitiAt - cot (3) m'^^ + + D*j^ ^ ' 



Dl = cos2/3M| [+^-^sm^ 9w^ , D^j^ = cos2/3Af| (^+'^sin^ew 



Dl = cos2/3M| ( + -sin^eiw 1 , D^^ = cos2/3Af| [--sin'ew]- (B5) 



APPENDIX C 



After defining the following product of propagators: 

D 



9 



(CI) 
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the integrals in Eqs. (jSJ and (|10|l can be cast in the form 



I til — 



Iti2 — 



(27r)4 [fc2- 



3" 



D 



fc2 — m: 

I ■ k 



dH f d'^k D 
d^k 

(2^ J (27r)4 [fc2-m2] \P-ml k^ ^ m: 



P — TO? 



for the case in which a photon/gluon is emitted by the t quark; 



Ihh = 



lHi2 — 



d'^k 



dH 



D 



(27r)4 J (27r)4 \k^ 



lA{l-kf~PP 



TO? 



fc2 



d'^k 



dH 



D 



I ■ k 



(2^)4 J (2^)4 [fc2-TO^±] |A:2-m2 ^2 _ ^2 
for the case of a photon emitted by the charged Higgs boson; 



d'^k 

d^k 

J2^ 



dH 
(2^ 



dH 

(2^ 



D 



{I + ky 



I kf - ;2fc2 + 3fc2(? . k) 4(Z • /c)2 - + ■ k) 

3 I F-TO? 



;2- 



l-{l + k) k-{l + k) 



(/ + A:)2 - 



;2 _ 



for the case of a photon/gluon emitted by the t squark; 



I HI 



d'^k 



dH 



D 



(2^)4 J (2^)4 [P-ml 



fc2 



14(rfc)2-?2fc2 
3 fc2_^2 



d4fc 



d4? 



Z2 



; • k 



(27r)4 7 (27r)4 [P - ml] P 



for the case of a photon/gluon emitted by the s squark; 



hi 



h2 



d^k 



dH 



D 



(27r)4 J (27r)4 



2ml{l-k) 14(rfc)2-;2fc2 
P - TO? ~ 3' 



/fc2 



d4fc 



d4? 



/ • A: 



2/2 



(27r)4 J (27r) 



P — 771? 



fc2 



P-TO? 



(C2) 
(C3) 

(C4) 
(C5) 

(C6) 
(C7) 



(C8) 
(C9) 



(CIO) 



(CU) 



for the case of a gluon emitted by the gluino g. 

These integrals were obtained from those corresponding to the diagrams (a)-(e) in Fig.|21after an expansion on the 
external momenta and a reduction of all tensorial structures to scalar ones. They can be further reduced to linear 
combinations of the two-loop vacuum integrals 



G{mi,pi;m2;p2;m3,p3) = 



d'^k 



dH 



1 



{2t:Y J {2nY [(; + A;)2 - to2]^^ [/2 - 7772]^^ [fc2 



-,2lP3 ' 



(C12) 



whose general solution is explicitly reported in Ref. |23l |. 

In the nondecoupling approximation of Refs. (iQUllj . only the integrals Iti2 and 1^12 are needed. They are evaluated 
after substituting {I + fc)2 with P in the expression for D. The terms proportional to / • fc in the curly brackets of 
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Eqs. 



Sj and ljC5|l are then dropped and the integrals are factorized into two one-loop integrals as 

d*k -2fc2 r dH 1 



^ti2 nondcc 



Ih 



i2 nondcc 



(2^)4 [fc2_^2]^ [p_^2^^] y (27r)4 



P — TIT- 


[P-mj] 


P-mj 









(2^)4 



;2 _ ^2 


[P-mj] 


P — m? 









(C13) 
(C14) 
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